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\  ABSTRACT 

Newton's  method  is  a  well  known  and  often  applied  technique  for 

computing  a  zero  of  a  nonlinear  function.  By  usinn  the  theory  of  generalized 

equations,  a  Newton  method  has  been  developed  to  solve  problems  arising 

in  both  mat  hem.it.i  a  1  programming  and  mathvmatii  al  economics^  ****** t*»»*rp!i¥  |H'H 
o-£jf  aA 

W.i  nlur  .mow  Newton's  method  for  generalized  equations  can  be  applied 

A, 

to  the  economic  equilibrium  problem  of  the  Project  Indej>endence  Evaluation 
System  (PIES)  Energy  Model.  Solutions  to  a  simplified  version  of  PIES 
arc  obtained  using  a  Newton  method,  and  comparisons  arc  made  to  solutions 
which  have  apjt'ured  in  the  1 1  terat  ure,  (Hogan  |9J). 
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k;nifk-ance  and  explanation 


Newton's  method  i  u  well  known  and  often  applied  technique  for 
computing  a  zero  of  a  nonlinear  function.  By  using  the  theory  of 
generalized  equations,  a  Newton  method  has  been  developed  to  solve 
problem:,  arising  in  both  mathematical  programming  and  mathematical 
economics  (Josephy  111)1. 

Recent  problems  of  practical  importance  in  the  study  and  computation 
ot  equilibria  for  large-scale  economic  systems  motivate  increased  interest 
in  the  use  of  Newton's  method  for  analyzing  and  computing  economic 
equilibrium.  -'no  such  problem  is  the  economic  equilibrium  of  the  Project 
Independence  Evaluation  System  (PIES)  Energy  Model.  Developed  by  the 
Federal  Energy  Administration  to  model  the  energy  sector  of  the  United 
States  economy,  PIES  formulates  an  equilibrium  problem  between  the 
suppliers  and  consumers  of  energy.  This  problem  can  be  modeled  as  a 
generalized  equation,  and  thus  can  be  analyzed  using  the  methods  developed 
in  Josephy  (11). 

We  present  compu tat lona 1  results  of  Newton’s  method  applied  to  the 
economic  equilibrium  problem  of  the  Project  Independence  Evaluation  System 
(PIES)  Energy  Model.  Solutions  to  a  simplified  version  of  PIES  are 
obtained  using  a  Newton  method,  and  comparisons  are  made  to  solutions 
which  have  appeared  in  the  literature. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
sunriary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


A  NEWTON  Mf.THOD  FOR  THE  PIES  ENERGY  MODEL 


Norman  H .  Josephy 


1 •  Introduction 

A  general ized  equation  is  a  set-valued  analogue  of  a  single-valued 
equation.  Many  problems  arising  in  mathematical  programming  and 
mathematical  economics  can  be  analyzed  by  generalizing,  to  the  set-valued 
ase,  classical  techniques  used  in  the  study  of  single-valued  equations. 
Preliminary  investigations  of  the  properties  and  applications  of 
generalized  equations  have  been  carried  out  by  Robinson  116-20). 

Josephy  (11,12)  develops  Newton  methods  for  computing  the  solution  of 
generalized  equations.  This  paper  presents  an  application  of  generalized 
equations  to  the  economic  equilibrium  problem  of  the  Project  Independence 
Evaluation  System  (PIES)  Energy  Model.  Developed  by  the  Federal  Energy 
Administration  to  model  the  energy  sector  of  the  United  States  economy, 
PIES  formulates  an  equilibrium  problem  between  the  suppliers  and 
consumers  of  energy.  This  problem  can  be  modeled  as  a  generalized 
equation,  and  thus  can  be  analyzed  using  the  methods  developed  in 
Josephy  (11,12). 

Section  2  begins  with  a  general  discussion  of  the  equilibrium 
problem  generated  by  trying  to  balance  supply  and  demand  of  consumable 
goods  in  a  competitive  market  environment .  An  example  published  by 
Hogan  (91  illustrating  the  major  aspects  of  the  PIES  model  is  detailed 
in  Section  3.  A  description  of  the  solution  procedure  used  by  the 
Federal  Energy  Administration  (FEA)  in  solving  the  equilibrium  problem 


Sponsored  by  the  Mni ted  Stetee  Army  tinder  Oontrect  No.  DAAG29-75-C-0024. 
This  material  is  based  upon  work  supported  by  the  National  Science  Founda¬ 
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chool  of  the  University  of  Wi sconsin-Madison. 


i  siven  in  Section  4.  Newton's  method  for  generalized  equations  is 
apj  lied  to  th<-  economic  equilibrium  problem  in  Section  5.  The  usi  ol 
a  Newton  method  to  compute  economic  equilibrium  is  first  suggested  in 
Robinson  1161.  An  alternative  approach  to  applying  a  Newton  method 
to  the  PIES  model  can  be  found  in  Eaves  [51 . 

Section  6  concludes  the  paper.  We  present  computational  results 
comparing  solutions  published  by  Hogan  and  the  solutions  obtained  by 
using  Lemke's  ulgoritlim  to  solve  the  linear  complementarity  subproblem:, 
of  the  Newton  iterations.  The  data  for  Hogan's  example  are  reproduced 
in  an  appendix. 

2.  Equilibrium  in  a  Market  Economy 

The  competitive  market  which  underlies  the  PIES  model  consists  of 
two  classes  of  agents,  suppliers  and  consumers,  and  two  classes  of 
goods,  factors  of  production  and  consumable  goods.  The  suppliers, 
faced  with  a  perceived  demand  for  consumable  goods,  convert  factors  of 
iroduction  into  consumable  goods,  sustaining  a  cost  for  conversion  to 
and  delivery  of  consumable  goods,  and  charging  a  price  for  those 
goods.  The  consumers  purchase  consumable  goods  at  levels  dependent 
upon  the  prices  of  all  goods. 

Suppose  there  are  n  consumable  goods.  bet  q  t  R^  and  p  c  R*', 
where  q^  is  the  consumption  level  of  the  ith  qood,  and  is  the  price 

of  the  ith  good.  For  a  perceived  demand  q  of  consumable  goods,  the 
supplier  will  charge  prices  no  lower  than  p  (q)  r  R°.  The  consumer, 
desiring  level  q  of  consumable  goods,  is  willing  to  pay  no  more  than 
prices  pD<q)  r  .  A  market  equilibrium  is  a  quantity  vector  q  and 
a  price  vector  p  such  that  p  ■  pD<q)  ■  pg(q).  That  is,  at  price  p. 


the  market  for  all  poods  will  clear.  The  producers  will  supply  at 
price  p  a  level  of  consumable  goods  q  which  the  consumers  are  willing 
to  purchase  at  price  p.  The  supply  relation  ps  is  typically  determined 
by  the  solution  of  a  cost  minimization  problem  modeling  the  engineering/ 
technological  processes  involved  in  conversion  of  factors  of  production 
to  delivered  consumable  goods.  The  demand  relation  pD  is  tradition¬ 
ally  a  behavioral  model  economet rical ly  determined  from  historical 
data.  The  following  example  will  illustrate  the  determination  of 
the  supply  relation  p<.. 

Example  1.  The  supplier  converts  two  factors  of  production 
into  two  consumable  goods.  The  factors  cost  dollars  per  unit, 
i»l,2,  respectively.  The  total  cost  of  converting  units  of  factor 
one  and  x,  units  of  factor  two  in  c^x^  ♦  c^x^.  *c  aSMunc  a  linear 
technology,  that  i^  for  factor  levels  Xj  and  x^,  an  amount 

a  1 1  x i  ♦  apx2  8°°^  onc  an^  an  •unount  a->ix]  *  a->->*2  8ooc*  lwo 

will  be  produced,  he  assume  there  is  a  demand  of  q  units  of  good 
onc  and  1  unit  of  good  two.  The  cost  minimization  problem  faced  by 
the  supplier  is  the  following  linear  program: 

minimize  CjXj  ♦  c^  subject  to 

allXl  4  a12X2  ”  q  • 

*21*1  4  *22*2  ‘  1  * 

Xj  •  0  ,  x2  ■  0  . 

The  dual  linear  program  is  given  by 


maximize  qPj  ♦  p.,  subject  to 

allPl  4  a21p2  "  C1  ' 

a  1 2P 1  4  a22P2  “  c2  * 

Pj  i  0  ,  P2  -  0  . 

The  dual  problem  is  illustrated  in  Figure  1.  It  is  assumed  that 
the  supplier  will  price  the  consumable  goods  at  marginal  price  levels 
Pj,  p^.  The  dual  linear  objective  function  is  a  line  qp^  ♦  p.,  in 
the  Pj  -  p^  plane.  For  q=0,  the  optimal  dual  solution  is  the  extreme 
point  1  ,  (O.p^).  This  remains  optimal  as  q  is  increased  until  the 
objective  function  is  parallel  to  face  1  -  2  .  At  that  value 

of  q,  every  point  on  the  line  segment  between  (0,p^)  and  (p^.p^)  is 
a  dual  optimal  solution.  As  q  is  increased  further,  extreme  point 


Figure  l.  Supplier's  Dual  Linear  Program 


2  ,  (Pj.i^)  becomes  the  unique  dual  optimal  solution.  This  behavior 
is  rcpeatcJ  with  face  2-5  when  q  is  increased  sufficiently  to 
make  the  objective  function  parallel  to  it.  For  larger  values  of 
q,  extreme  point  3  ,  (Pj,0)  remains  the  unique  dual  optimal  solu¬ 

tion.  The  relationship  between  q  and  p^  is  illustrated  in  Figure  2. 
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Figure  2.  Supplier's  Price  vs.  Quantity  Curve  pg 

The  consumer's  demand  curve  p^  is  obtained  by  fitting  historical 
data  to  prcspccificd  functions  (see  FEA  (61,  Halvorscn  18]).  In 
particular,  a  linear  relationship  between  the  logarithm  tn(q^)  and 
the  logarithms  tn(p^)  is  postulated.  Thus,  for  the  case  of  n  consum¬ 
able  goods,  the  demand  relation  is  determined  by  the  following  set 
of  equations: 

»n(qj/q^)  -  CjjtnCPj/p®)  ♦  •••  ♦  «ln  *n(Pn/P°) 

tn (q  /q°)  ■  c  .  ln(p,/p?)  ♦  ♦  e  ln(p  /p°)  , 

v^n  ^n'  nl  'rrM  nn  vrn  *ir 

where  c  is  the  elasticity  matrix  and  q^,  p*^  arc  fixed  reference  values 

of  pD,  that  is,  PD(q°)  -  p°  . 

The  market  equilibrium  for  n=2  is  illustrated  in  Figure  3. 

It  should  be  noted  that  p*,  the  equilibrium  price  vector,  is  not  at  a 
value  given  by  an  extreme  point  of  the  supplier's  dual  linear  program. 
Thus,  if  q  were  assigned  the  value  q*  in  the  supplier's  cost  minimiza¬ 
tion  problem  and  the  simplex  method  was  used  to  solve  the  resulting 


linear  program,  one  would  obtain  either  p  or  p  as  optimal  dual  variable^ 
and  not  the  equilibrium  value  p*. 


q* 

Figure  3.  Market  Equilibrium 

For  the  case  of  n  consumable  goods,  the  equilibrium  quantity 

q*  and  price  p*  and  factor  levels  x*  for  the  competitive  market  with 

supply  modeled  as  a  linear  program 

minimize  <c,x>  total  cost  of  production 

Ax  *  q  demand  requirements 

Bx  *  b  non-demand  constraints 

x  "  0  non-negative  factor  levels 

and  a  log- linear  consumer  demand  model 

*n(q  ■  £  c  ln(p  /p°)  ,  i*l , . . .  ,n 

11  j-1  J 

is  a  triple  (q*,p*,x#)  i  0  satisfying  the  following  equilibrium  condi¬ 
tions. 


(EQ.  1)  x*  solves:  minimize  <c,x>  subject  to 

Ax  «  q*  ,  Bx  =  b  ,  x  =  0  . 

(EQ.  2)  For  some  multiplier  s*,  (p*,s*)  solve  the  dual 

problem:  maximize  <q*,p>  ♦  <b,s>  subject  to 

T  T  <  > 

A  p  ♦  B  s  «  c  ,  (p,s)  *  0  . 
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(IQ.  3) 


mCq-A^) 


n 

Z 


j  =  l 


e.,  tn(p*/n  ) 

ij  *Y 


i-1 , . . . ,n  ■ 


In  non-technical  terms,  the  equilibrium  conditions  put  a  constraint, 
the  demand  relation,  on  the  optimal  dual  solution.  The  constraint 
relates  the  optimal  dual  solution  to  the  cost  coefficients  of  the 

dual  linear  program.  He  note  that  it  is  incorrect  to  solve  the 

equilibrium  problem  by  including  the  demand  constraints  in  the  sup¬ 
plier's  linear  program.  The  following  example  illustrates  this  point. 

Example  2.  Consider  the  supply  linear  program  minimize  x 
subject  to  x  *  q,  x  =  0  ,  and  its  associated  dual  program 

max  p  q  subject  to  p  »  1  ,  p  *  0  . 

The  optimal  dual  price  is  p*  =  1  for  any  q  >  0  .  Suppose  the  demand 

relation  is  given  by 

4(1 -q)  if  0  *  q  “  1 

PpCq)  -  i 

0  if  q  >  1 


The  equilibrium,  as  shown  in  Figure  4  is  p*  ■  1,  q*  *  1. 


7 


Figure  5.  The  Hogan  PIES  Example. 


1  unit  of  capital  and  1  ton  of  steel  per  ton  of  coal,  and  will  cost 
S  dollars  per  ton  of  coal.  The  variable  c^  represents  this  level 
of  coal  production.  Additional  production  above  300  tons  up  to 
600  tons  of  coal  consur.cs  5  units  of  capital  and  2  tons  of  steel  per 
ton  ot  coal,  and  costs  6  dollars  per  ton.  This  level  of  coal  produc¬ 
tion  above  300  tons  is  denoted  by  Cj,.  Between  600  and  10C0  tons 
of  coal  mined,  each  tor.  cf  coal  requires  10  units  of  capital,  3  units 
of  steel,  and  costs  8  dollars  per  ton.  This  level  of  production  is 
denoted  by  c^.  Thus,  the  total  tonnage  of  coal  mined  in  region  1 
is  ♦  Cp  ♦  c^j  ,  where  0  •  c  »  300  ,  0  ■  c^  ■  300  • 

0  ■  Cjj  ■  400.  The  total  cost  of  this  production  is  5c ^ ^  ♦  bc^.,  ♦  8c 

and  requires  1  *  c  ^  ^  ♦  2*c^  ♦  3*Cj^  tons  of  steel  and 

1‘Cjj  ♦  S'Cj^  *  *®*C13  capital.  The  cost  curve  for  coal 

■ined  in  region  1  is  given  in  Figure  1  .  A  similar  analysis  applies 
to  the  other  coal  and  oil  regions.  Data  for  those  regions  arc  given 

in  the  .i;  ;  »>nd  l  x  . 


Figure  *  . 


Coal  if  transported  at  level  cT. .  from  coal  region  i  to  demand 


region  j  at  a  cost  per  ton  given  in  the  appendix.  Crude  oil  is  trans 

ported  from  oil  region  i  to  refinery  j  in  a  quantity  oT^ ,  where  it 

is  refined  into  a  fixed  proportion  of  light  and  heavy  oil.  Refinery 


1  converts  60\  of  its  crude  oil  into  light  oil,  while  refinery  2  con 


verts  50*  of  its  crude  oil  into  light  oil.  Light  oil  L. .  and  heavy 


oil  H  arc  transported  from  refinery  i  to  demand  region  j.  The 


transportation  costs  for  crude  oil  oT  and  refined  oil  L,  H,  are  given 


in  the  appendix 


The  linear  program  representing  the  minimum  cost  of  production 
and  delivery  of  consumable  goods  to  meet  a  demand  q  c  R*’  is  of  the 


demand  requirements 


transportation  flow  balance 
factor  of  production  upper  bounds 


steel  aqd  capital  resource  constraint 
non-negitive  factors  of  production 


The  components  of  x  represent  the  levels  of  production  of  coal  and 
crude  oil,  levels  of  coal  transported  to  demand  regions,  levels  of 
cruJe  oil  transported  to  refineries,  and  levels  of  refined  oil  trans 


ported  to  the  demand  regions.  The  coefficient  r.  is  the  unit  cost 


of  the  activity  associated  with  x  .  The  vector  Ax  consists  of  levels 


of  consumable  goods  delivered  to  the  demand  regions.  Bx  represents 


the  conservation  of  materials  transported  through  a  region  or  re¬ 
finery.  Fx  is  the  vector  of  steel  and  capital  consumed  in  the 
processing  of  factors  x.  The  equations  and  inequalities  given  by 


this  linear  program  are  detailed  in  the  appendix. 

The  demand  relationship  between  the  desired  quantity  of  a  con¬ 
sumable  good  in  a  particular  demand  region,  q^,  and  the  prices  con¬ 
sumers  are  willing  to  pay  for  such  goods,  {p^}  j 6,  is  given 
by 


^(q^) 


6 

Z 

j-1 


e.  tn(p  /p  ) , 
ij  J  V 


i- 1 ....  ,6 


The  values  of  the  data  q 


0  j 

p  and 


arc  given  in  the  appendix. 


4  .  KKA  So  Hit  Ion  A I  gor  1 1  hr 

The  motivation  for  the  FEA's  solution  algorithm  for  computing 
the  market  equilibrium  in  the  PIES  model  (Hogan  (9),  FF.A  (7)  )  fol¬ 
lows  from  the  observation  that  finding  a  :ero  of  a  gradient  mapping 
can,  under  suitable  conditions,  be  replaced  by  minimizing  a  scalar- 
valued  potential  function  defined  in  terras  of  the  line  integral  of 
the  gradient  mapping  (Ortega  and  Rhcinboldt  (if)  ).  This  observation 
can  be  applied  to  the  market  equilibrium  problem  in  the  following 
manner.  The  kuhn-Tucker  conditions  for  the  suppliers  linear  program 
(P)  with  perceived  demand  q*, 


(P)  min  <c,x>  subject  to  Ax  «  q*  ,  Bx  *  b  ,  x  ■  0  . 

are  the  classical  complementarity  conditions  (Luenbergcr  |14]  ): 
there  exists  x*.  p*.  s*,  such  that 
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The  demand  constraint  is  given  by 


Note  that  the  equilibrium  conditions  (EQ.  1)  and  (EQ.  2)  of  the  prev 


ious  section  are  equivalent  to  conditions  (1),  (2),  and  (3),  while 


equilibrusn  condition  (EQ.  3)  is  equilvalcnt  to  equation  (4).  Thus 


the  four  conditions  (1)  -  (4)  determine  equilibrium  values  of  con 


sumption  quantity  q* ,  price  p#,  production  factor  levels  x 


Suppose  p.  is  a  gradient  mapping 


demand  constraint  multipliers  s 


that  is,  for  some  function  g:R  ♦  R 


The  Kuhn-Tucker  conditions  for  the  non-linear  minimization  problem 


minimize  <c,x>  -  g(q)  subject  to 


are  the  following 


(3') 


Bx*  -  b  *  0  ,  s*  =  0  ,  <s*,  Bx*-b>  *  0 
(4')  -Vg(q*)  ♦  p*  =  0 


Using  (5)  to  substitute  p()(q*)  for  Tg(q*)  in  (4'),  we  see  that  (l1), 
(2'),  (3'),  and  (4')  arc  identical  to  (1),  (2),  (X),  and  (4).  Thus, 
under  the  assumption  that  p()  is  a  gradient  mapping,  we  can  replace 
the  equilibrium  problem  with  finding  a  Kuhn-Tucker  point  for  the  non 
linear  programming  problem  (NLP).  However,  p;)  is  a  gradient  mapping 
in  a  neighborhood  of  a  point  q*  if  and  only  if  the  derivative  (q) 
is  symmetric  for  q  in  that  neighborhood  (Ortega  and  Rheinboldt  1151, 
p.  95).  This  condition  does  not  hold  for  the  P1LS  model  p()  function. 
To  circumvent  this  difficulty,  the  algorithm  used  in  FLA  (17)) 
employs  a  zero  cross-elasticity  approximation.  That  is,  the  original 
demand  relation 

0  n  0 

(6)  tnq  -  tnqt  ♦  l  e  tn  (p  /p  )  i«l . n 

1  1  j-1  **  J  J 


is  replaced  by 


(«■') 


tnqj  ■  tnq°  ♦  c..  tnfp^pj) 


i*l , 


,n 


This  yields  the  consumer  demand  function  p^  with  components 

0  0  Ki* 

(7)  PD(<lli  *  pi  iBl . n 

which  has  a  symmetric  derivative,  h'c  are  now  faced  with  solving 
problem  (NLP),  with  g  given  by 
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n  arc  negative 


monotone  non-increasing,  and  we  can  thus  approximate  (-g)  by  a  piece 


wise  linear,  convex  function  whose  value  at  any  q  can  be  represented 


ns  the  solution  of  a  linear  program  with  q  appearing  as  a  constant 


in  the  constraints.  Thus,  replacing  (-g)  in  (NLP)  by  a  linear  pro 


gram,  we  arrive  at  a  large  linear  program  of  the  form 


which  is  an  approximation  to  (NLP),  whose  solution  yields  an  equi¬ 
librium  point  of  the  original  market  equilibrium  problem.  The  details 
of  the  approximations  and  formulations  discussed  above  can  be  found 


in  Hogan  (9]  ,  and  FEA  171 


The  FEA  solution  algorithm  proceeds  as  follows.  The  linear 


programming  problem  (LP)  is  solved  to  obtain  an  estimate  of  the 


equilibrium  price  vector  p 


is  computed  by  setting  p  >  p 


elasticity  approximation  is  made  about  (q 


placed  by 


1  .  Newton' s  Met  hod  for  TIES  Equilibrium 

The  market  equilibrium  for  the  PIES  model  is  a  quadruple 
(x*,p*,q*,s*)  satisfying  conditions  (1),  (2),  (3),  and  (4).  Let 
F:  Rl*n  -*  Rl*n  be  defined  by 


(9)  F(x,p,s,q):  * 


1 

c 

0 

-at 

-bt 

-s 

0 

X 

o 

♦ 

A 

0 

0 

-I 

P 

-b 

B 

0 

0 

0 

s 

'PD(q). 

.  o 

I 

0 

0. 

.  q  . 

and  K:  «  r|  <  Rn  ,  where  t  «  m*n«v  ,  x  c  R1**  ,  p  c  Rn,  s  c  r'  ,  and 
q  c  Rn  .  Then  the  equilibrium  conditions  can  be  written,  using  the 
notation  of  Josephy  1111,  as  x*,p#,s*,q*  solving 


(10)  0  c  F(x,p,s,q)  ♦  Nj, (x ,p, s ,q) 

The  linearization  of  (10)  about  x,p,s,q  is 


(11) 


0  c 


c 

,  T  T  -i 

0  -A1  -B*  0 

X 

0 

AO  0  -I 

P 

♦ 

. 

-b 

BO  0  0 

s 

t  . 

•pD(qi*pD(q)q 

0  1  0  -p£(q)_ 

q 

♦Nk(x,p,s,q) 


This  corresponds  to  the  linear  complementarity  problem  obtained  from 
conditions  (1),  (2),  (3),  and  (4)  by  replacing  (4)  with 
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6.  Computational  Results 


It  has  been  shown  in  Josephy  (11)  that  Newton's  method  for  generalised 


equations,  when  applied  to  the  nonlinear  complementarity  problem 


results  in  a  sequence  of  linear  complementarity  problems.  We  have 


shown  how  the  economic  equilibrium  problem  of  the  PIES  energy  model 


r>  iits  in  a  nonlinear  complementarity  problem.  Thus,  to  compute 


equilibrium  solutions  for  PIES,  we  solve  linear  complementarity  problems 


The  solution  technique  we  will  use  is  Lemke's  complementary  pivot 


algorithm  (Cottle  and  Dantzig  (4)).  Not  all  linear  complementarity 


problems  are  solvable  by  Lemke's  algorithm.  However,  it  is-  shown  in 


Josephy  [131  that  for  the  data  given  by  Hogan,  the  linear  complementarity 


problems  resulting  from  applying  Newton's  method  for  generalized  equa 


tions  to  Hogan's  example  are  solvable  by  Lemke's  algorithm 


Table  1  summarizes  the  equilibrium  quantities  and  prices  for 


Hogan's  example.  Column  1  contains  Hogan's  published  equilibrium 


values.  Column  2  contains  the  equilibrium  values  computed  by  Newton's 


method  as  developed  here,  with  Lemke's  algorithm  applied  to  the  linear 


complementarity  problem.  Hie  equilibrium  values  in  Column  3  are 


computed  as  were  those  of  Column  2,  except  that  the  off-diagonal 


elements  of  the  elasticity  matrix  are  set  to  zero.  In  each  case 


in  which  Newton's  method  was  applied  to  the  equilibrium  nonlinear 


complementarity  problem,  convergence  occurred  in  three  iterations 


TABLE  1 

EQUILIBRIUM  QUANTITIES  AND  PRICES 


Hogan 

Newton 

Newton  with  Zero 
Cross-Elasticities 

ql 

996 

1016 

1029 

^2 

910 

912 

913 

^3 

120S 

1201 

1212 

% 

1229 

1223 

120S 

<»S 

996 

998 

1012 

% 

1020 

1011 

992 

Pi 

11.3 

11.7 

11.6 

P2 

13.4 

13.7 

13.6 

Pj 

IS. 4 

IS. 8 

15.7 

P4 

IS. 6 

16.0 

IS. 9 

Ps 

11. S 

11.9 

11.7 

Pft 

12.0 

12.4 

12.7 
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with  <in  error  less  than  10  The  difference  in  the  values  between 
column  l  and  column  2  are  due  to  the  error  introduced  in  the  FEA 
algorithm  when  a  linear  programming  approximation  (LP)  is  made  to  the 
nonlinear  programming  problem  (NLP).  Ahn  [1]  has  shown  that,  under 
suitable  conditions,  the  FEA  algorithm  will  converge  to  an  equilibrium 
if  (NLP)  rather  than  (LP)  is  solved  at  each  iteration. 

Table  2  lists  the  equilibrium  values  of  the  problem  variables. 

The  notation  used  in  both  tables  is  defined  in  the  appendix. 

Rows  1  and  2  of  Table  3  list  the  optimal  dual  prices  for  the 
supplier's  linear  program  when  the  fixed  demand  q  is  set  to  the  Hogan 
equilibrium  and  to  the  Newton-computed  equilibrium.  The  third  row 
of  Table  3  lists  the  dual  prices  for  the  supplier's  linear  program 
with  fixed  demand  q  set  to  the  Newton  equilibrium  with  q^  decreased 

by  .1  percent.  These  numbers  were  obtained  by  applying  the  simplex 
method  to  the  supplier's  linear  program. 

The  numbers  in  Table  3  indicate  that  the  dual  feasible  region 
has  extreme  points  whose  price  components,  as  listed  in  rows  1  and 
2  of  the  Table,  are  not  close.  However,  a  small  change  in  only  one 
component  of  the  dual  program  objective  coefficient  vector  causes 
the  optimal  prices  to  jump  between  the  two  extreme  points.  Also 
note  that  neither  dual  solution  Is  the  equilibrium  price  vector. 

This  type  of  occurrence  is  discussed  in  Section  2. 

We  remark  in  closing  that  the  generalized  equation  framework 
is  eminently  capable  of  modeling  equilibrium  problems  containing 
equations  or  unconstrained  variables.  An  example  of  a  more  general 
equilibrium  problem  than  that  discussed  in  this  naper  is  the  following. 
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1  1 1,111  t  ■'  ^^11 


n 


TABLE  3 


DUAL  LINEAR  PROGRAM  PRICES 


Pi 

P2 

P3 

P4 

PS 

P6 

Hogan 

8.5 

10.5 

16.4 

16.6 

2.6 

3.1 

Newton 

37.6 

39.6 

45.4 

45.6 

41.6 

42.1 

Newton  with 
.  1\  Change 

8.5 

10.5 

16.4 

16.6 

2.6 

3.1 

It  is  of  interest  to  gauge  the  effect  of  the  energy  sector  of  the 
United  States  economy  on  the  economy  as  a  whole  (see  Ask  in  12,31  and 
Hogan  and  Manne  (10)).  One  apj  roach  to  measuring  this  effect  is  to 
include,  in  the  microeconomic  energy  model,  macroeconomic  equations 
relating  the  level  of  consumer  purchases  both  to  prices  of  goods  and 
to  macroeconomic  variables  such  as  gross  national  product.  These 
additional  equations  can  be  i ncorporated  into  the  equilibrium  formula¬ 
tion  and  the  techniques  discussed  in  Josephy  111)  can  be  applied  to  the 
generalized  equation  formulation  of  the  enlarged  equilibrium  problem. 
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DATA  FOR  HOGAN'S  EXAMPLE 


coal  in  coal  region  i  at  increment  j 


oil  in  oil  region  i  at  increment  j 


coal  transported  from  coal  region  i  to  demand  region  j 


oil  transported  from  oil  region  i  to  refinery  j 


light  oil  transported  from  refinery  i  to  demand  region  j 


heavy  oil  transported  from  refinery  i  to  demand  region  j 


TABLE  1 


QUANTITY  AND  PRICE  NOTATION 


Heavy  Oil 


Heavy  Oil 


The  following  data  specify  the  supply  side  of  the  economic 


equilibrium  problem 


Coal  production  level  limits: 


coal  region  1:  0  <  Cjj  <  300  ,  0  <_  f.  300  ,  0  £  Cjj  <_  400 

coal  region  2:  0  <  c7j  <  300  ,  0  _<  c^j  *  300  ,  0  _<  c^j  £  t>00 

Oil  production  level  1 i m its: 

oil  region  1:  0  0 ^  ^  <  1100  ,  0  0 ^ ^  2L  1200  ; 

oil  region  2:  0  £  0,^  *c  1300  ,  0  <  0^  j<  1100  . 

T  r  an sport  at  ion  f low  balance: 


coal 

regions: 

cil  *  ci2  •  ci3  -  cTil  •  cTi2 

-  cTi3 

oil 

regions : 

°ll  *  °12  '  °TU  ‘  °Ti2  ’  ° 

i-1.2 

Refinery  yields: 

Refinery  1  Refinery  2 

Light  Oil  .6  .5 

Heavy  Oil  .4  .5 

T ransportation  flow  balance: 


light  oil ,  refinery 

1: 

L11 

*  L12  ‘ 

.6 

f°Tll 

♦  °T21> 

-  0 

light  oil,  refinery 

2: 

L21 

^  L22  • 

.5 

C°T  2 

4  °T22> 

«  0 

heavy  oil,  refinery 

1: 

■u 

H12  * 

.4 

(°T11 

♦  °T21> 

•  0 

heavy  oil ,  refinery 

2: 

H21 

«22 

.5 

C°T12 

4  oT22> 

-  0 

Demand  Constraints: 


coal:  cTjj 

*  cI2!  i  ’l 

CT12 

♦  CT„  I  02 

light  oil: 

L11  *  L21  i  ’j 

lJ2  *  l22  i  ’4 

heavy  oil: 

"ll  *  H21  i  ’5 

H12  *  H22  i  \ 
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